Island nucleation in thin-film epitaxy: A first-principles investigation 
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We describe a theoretical study of the role of adsorbate interactions in island nucleation and 
growth, using Ag/Pt{lll) heteroepitaxy as an example. From density-functional theory, we obtain 
the substrate-mediated Ag adatom pair interaction and we find that, past the short range, a repulsive 
ring is formed about the adatoms. The magnitude of the repulsion is comparable to the diffusion 
barrier. In kinetic Monte Carlo simulations, we find that the repulsive interactions lead to island 
densities over an order of magnitude larger than those predicted by nucleation theory and thus 
identify a severe limitation of its applicability. Copyright 2000 by The American Physical Society. 
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Island nucleation is often the first step in thin-film epi- 
taxy and is, thus, relevant to the synthesis of a wide 
variety of interfacial materials. Achieving a quantita- 
tive understanding of the island morphologies (i.e., sizes, 
shapes, density, spatial distribution, etc.) that develop in 
the initial stages of thin-film growth is also important for 
fundamental reasons. Since thin-film epitaxy frequently 
occurs away from equilibrium, the kinetics of deposition 
and surface diffusion play a key role in governing island 
morphology and there is great variety in the resulting 
structures. Considering shapes for example, islands 
can be fractal-like or compact and triangular, hexagonal, 
square, rectangular, etc.. Each of these structures is a 
signature of an intricate kinetic balance and reflects a 
complex set of interatomic interactions that is unique for 
each material. 

Despite the complexity and potentially enormous vari- 
ety in growth morphologies, certain aspects of island nu- 
cleation and growth appear to be common to many differ- 
ent systems. In a general description, gas-phase species 
are deposited onto an initially bare solid substrate with 
a rate F. These species hop on the surface with a rate 
D = i'Qe~^f'/^^'^ , where vq is the preexponential factor, 
E^^ is the diffusion-energy barrier for an isolated species, 
/cs is Boltzmann's constant, and T is temperature. Hop- 
ping mediates the aggregation of adspecies into nuclei, 
which either dissociate with an energy barrier Ed,i^ if 
they are below a critical size i, or grow subsequently to 
become stable islands. Initially, the formation of island 
nuclei is the main process taking place. As the surface 
coverage increases, it becomes increasingly likely that de- 
posited species will add to stable islands and promote 
their growth instead of forming new nuclei. These gen- 
eral features can be captured in a mean-field theory for 
the stable island density ||]. In the island growth 
regime, this expression has the form 



(1) 



Although the utility of a general expression cannot be 
overstated, Eq. ^ cannot describe all aspects of thin-film 
epitaxy and it is important to understand its limitations. 
In the interest of achieving a complete and predictive 
model for thin-film morphology, it is clearly desirable to 
have an approach that is as free as possible from arbi- 



trary parameters or assumptions. In this Letter, with 
an aim toward this ideal approach, we present the re- 
sults of a combined kinetic Monte Carlo (kMC) and first- 
principles, density-functional theory (DFT) study of is- 
land nucleation in a model for the growth of Ag on a 
monolayer (ML) of Ag on Pt(lll). Our choice of this 
model system was motivated by intriguing results from 
recent, low-temperature, scanning-tunneling microscopy 
(STM) studies 0-0, in which Eq. [l]was used to obtain 
the energy barrier and preexponential factor for adatom 
hopping. Shown in Table I are the parameters obtained 
in these studies. 

A striking feature of the experimental results is that 
the preexponential factors are significantly smaller than 
would be anticipated for systems such as these. For ex- 
ample, from ah initio calculations, Ratsch and Scheffler 
p3[ find a preexponential factor oi vq — l.i x 10^^ s~^ 
for a Ag adatom on 1-ML-Ag/Pt(lll), with a diffusion 
barrier of = 63 meV. Inserting the experimental and 
theoretical values for the diffusion parameters into Eq. |l| 
in the low-temperature limit where i = 1 and Ed.i = 0, 
we see that the experimental island densities are about 
an order of magnitude higher than predictions based on 
the theoretical diffusion parameters. Here, we investigate 
the origins of this discrepancy. Our DFT-kMC model in- 
cludes many features of the complex potential-energy sur- 
face experienced by Ag adatoms during thin-film growth 
and is free from several of the assumptions in Eq. 
We find that one of these assumptions - that interac- 
tions between adsorbed species do not extend beyond a 
short range - is violated. For systems with low diffusion- 
energy barriers [such as Ag on 1-ML-Ag/Pt(lll)], we 
show that these long-range, adatom-adatom interactions 
play an important and previously underestimated role in 
island nucleation and growth. 



System 


El (meV) 


(s-^) 


Ref. 
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Al on Al(lll) 


42 


8 X 10'' 




Ag on 1-ML-Ag/Pt(lll) 


60 


10" 





TABLE I. Experimentally determined diffusion-energy 
barriers and preexponential factors. 



1 



The DFT calculations [|4| are performed using the 
plane- wave, pseudopotential |15| method within the gen- 
eralized gradient approximation [ p^ . Previously, Ratsch 
et al. [0 showed in DFT calculations that the diffusion- 
energy barrier of an Ag atom on the 1-ML-Ag/Pt(lll) 
substrate is essentially the same as that on a strained 
Ag(lll) substrate, in which Ag is given the lattice con- 
stant of Ft. Thus, to model the heteroepitaxial system 
we use strained Ag(lll), in which the lattice constant is 
set to a value of 4.01 A. This value is 4.61 percent smaller 
than our calculated lattice constant for bulk Ag. We use 
the supercell approach to describe the surface, which is 
modeled as a (4 x 4 x 4) slab with a vacuum spacing 
of five interlayer distances. The cut-off energy is 50 Ry 
and we use 4 k points to sample the full surface Brillouin 
zone. The top layer of a bare slab is fully relaxed. Subse- 
quently, an adatom is placed in a binding site (fee and hep 
three- fold hollow sites), and its height is optimized with 
respect to the fixed substrate. To calculate adatom inter- 
action energies, two (or more) adatoms are placed on the 
relaxed (and fixed) substrate with heights fixed to values 
from the single-adatom calculations. In this way, we seek 
to isolate the electronic interaction between adatoms in 
binding sites. With simultaneous relaxation of both the 
adatoms and surface atoms, we can resolve the role of 
substrate-mediated, elastic interactions in the total in- 
teraction energy. Full relaxation of a few trial structures 
and inspection of the forces in our partially relaxed slabs 
indicates that elastic interactions are not highly depen- 
dent on adsorbate configuration and that our results will 
change by 10 meV or less with full relaxation. 

The total interaction energy AE for a periodic slab 
containing N adatoms, of which M are at binding site 
a and {N — M) are at binding site 6, is given by AE = 
E^S^Ij^ - ME-s^, -{N- M)El^^ + {N- l)Es. Here, 
^s'+N total energy of a slab with N adatoms, Eg^^ 

and Eg^^ are the total energies of slabs containing one 
adatom, and Es is the total energy of a bare slab. In 
the DFT supercell approach, the total interaction energy 
is comprised of interactions between different adatoms in 
the slab and interactions between adatoms in the slab 
and the periodic- image adatoms. We can express AE 
as a function of these interactions using the lattice-gas 
Hamiltonian approach (see, e.g. uM), which yields 



(2) 



Here, the summations run over all sites i in the slab and 
all sites j and k in the supercell (which includes both 
the slab and its periodic images), n„j is unity if site m 
(to = k) is occupied and zero, otherwise, F^^^(Rij) is 
the pair interaction between two adatoms on sites i and j, 
and y'^'^^(Rij , Ri_fc) is the trio interaction between three 
adatoms on sites i, j, and k. We neglect higher-order in- 
teractions. Another assumption implicit in Eq. H is that 



the interaction between adatoms at a fixed distance is 
independent of whether these atoms occupy fee or hep 
sites. We confirmed this assumption in one trial calcu- 
lation. Finally, the adatom binding energies on fee and 
hep sites are virtually equal: The fee site is favored by 
less than 3 meV. 

Thus, for a given adatom configuration, we express AE 
as a sum of pair and trio interactions with unknown co- 
efficients. From 18 different configurations, we obtain 
a system of linear equations and solve these for pair- 
interaction coefficients up to the 13*'^-neighbor, as well as 
for 5 different trio interactions. We assume that all other 
interaction coefficients are zero. To verify our parameter- 
ization of Eq. (2) , we used our interaction parameters to 
predict the total interaction energy in several, additional 
test structures. All of the predicted values agreed well 
with values from DFT calculations. 

The pair interaction is shown in Fig. ^ where we also 
show results for Ag on unstrained Ag(lll). For both 
surfaces, this interaction is strongly attractive at the 
nearest-neighbor distance and repulsive at longer dis- 
tances. It is interesting to consider the origins of the 
long-range repulsion. At these distances, the interaction 
could be due to substrate-mediated elastic interactions 
or of electronic origin jl^]. Since we find that elastic 
interactions play a small role here, the repulsion is pri- 
marily an electronic effect. Each adatom induces a small 
perturbation in the electron density, which decays with 
distance from the adatom in an oscillatory manner. The 
asymptotic tail, which is expected to decay with distance 
d as d~^ (or as d~^ , if a partially filled surface state is 
involved), is a Friedel-type oscillation. Friedel oscilla- 
tions have been imaged as concentric, ring-like, features 
around defects in low-temperature STM studies of several 
noble-metal surfaces [including Ag(lll)] ||2^. We expect 
the Friedel tail to extend to much longer distances than 
can be probed in DFT calculations. However, interac- 
tions associated with the Friedel tail should be weaker 
than those probed here. Thus, the central-ring interac- 
tion resolved here will have the most significant ramifi- 
cations for thin-film morphology. 

From Fig. 0(a), we see that the magnitude of the re- 
pulsive ring for Ag on strained Ag(lll) is comparable 
to the diffusion-energy barrier for an isolated adatom. 
For Ag on Ag(lll), it appears that the repulsive interac- 
tion is weaker and the diffusion-energy barrier is larger. 
The diffusion barriers reported here are obtained with 
full relaxation of both the first-layer slab atoms and the 
adatom. Our barriers are in good agreement with experi- 
mental values |l2) for Ag on 1-ML-Ag/Pt(lll) (60 meV) 
and on Ag(lll) (97 meV) and with those of Ratsch and 
Scheffler |l^. If the interaction energy and the diffusion 
barrier are of comparable size, we expect interatomic in- 
teractions to significantly inffuence adatom diffusion and 
island formation. Since Eq. ^ neglects the inffuence of 
long-range interactions, it is unclear if this expression is 
accurate under these circumstances. 

To resolve the effect of long-range interactions on thin- 
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FIG. 1. Pair interaction energy as a function of distance 
from a central adatom, shown in black, for (a) Ag on strained 
Ag(lll) and (b) Ag on Ag(lll). Also noted are the diffu- 
sion-energy barriers, E'j^, of isolated adatoms on these sub- 
strates. 



film growth, we developed a kMC model employing the 
general method of Fichthorn and Weinberg and in- 
corporating the pair potential for Ag on strained Ag(lll) 
shown in Fig. |](a). In the initial stages of thin-film epi- 
taxy, the surface coverage is low and pair interactions 
are likely to be the only significant interactions govern- 
ing island nucleation and growth ||2l[] . In our kMC model, 
atoms are deposited onto a fcc(lll) substrate with a rate 
of = 0.1 ML/s. An adatom hops from site i to site j 
with a rate given by Di—,j = i'Qe~^^^^^''^'^ , where Ei^j 
is the energy barrier to hop from site i to j. For the 
hopping-rate parameters, we use vq = 10^^ s^^ The 
energy barrier is given by Ei^j = E^ ■ — Ei, where Ei is 

the energy with an atom at site i and Ef ■ is the energy 
of the transition state between sites i and j. In general, 
eI ■ should depend on both Ei and Ej. Considering pos- 
sible permutations of adatom configurations with IS**^- 
neighbor interactions, ^ 10^"* different, diffusion-energy 
barriers could occur. To make the problem tractable, we 



adopt a simple model, in which Ei. 



El 



\{Ej-E, 



All of the quantities in this equation are obtained from 
DFT calculations. We have tested this equation for a trial 
geometry in which an adatom with four fee 9**^ neighbors 
hops to a nearest-neighbor hep site where it has two 7*'' 
and two 12'^ neighbors. From our simple model, we find 
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FIG. 2. Arrhenius plot of the island density as a function 
of temperature from an impermeable repulsive ring (squares) , 
the DFT-kMC model (circles), repulsive ring with e_B, = 25 
meV (diamonds), and nucleation theory (triangles). 



Ei^j = 53 meV, which is in remarkable agreement with 
the value of 46 meV from DFT calculations. 

We simulated thin-film epitaxy over temperatures 
ranging from 40-70 K and determined island densities 
in the beginning of the island growth regime. These low 
temperatures are in the range of the experimental studies 
(cf.. Table I). At such low temperatures, Eq. |l| reduces 
to the form ^ (F/D)^^'^. Fig. ^ shows an Arrhenius 
plot of the island density from our DFT-kMC model as 
a function of temperature. Also shown in Fig. ^ is the is- 
land density predicted by nucleation theory for the values 
of F, i/Q, and E^ used here. To quantitatively compare 
nucleation theory with the simulations, a proportionality 
coefficient rj is needed in Eq. {i.e., = r^{F/Df/^). 
This coefficient is related to the efficiency of the islands 
in capturing adatoms. Using a self-consistent approach, 
rj = 0.25 and values of -q ranging from 0.2 to 0.23 have 
been found in kMC simulations of Ag island nucleation 
on Pt(lll) H]. Here, we use -q = 0.25. 

In Fig. 2, we see that the DFT-kMC island densities 
are an order of magnitude (or more) above the theoreti- 
cal values. To understand this, we construct a caricature 
model, in which we replace the set of pair interactions 
shown in Fig. 0(a) with a nearest-neighbor attractive in- 
teraction and a uniform, repulsive ring of strength er 
at distances 10-13. By varying the magnitude of sr, we 
span the entire range of possible behaviors in this system. 
As cxD, the island density assumes a constant, max- 

imum value that is independent of temperature (cf.. Fig. 
^ . This is because island nucleation can only occur when 
one atom is deposited within the repulsive ring of another 
and it is governed by the temperature-independent depo- 
sition rate. In this regime, many adatoms are isolated by 
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repulsion in the initial stages of deposition. Each isolated 
adatom becomes a stable island when another atom is 
deposited into its ring and the resulting island density is 
significantly higher than in the absence of such a ring. 

As er is decreased, diffusing adatoms are increasingly 
able to surmount the ring barrier and a second channel 
for island nucleation and growth (via long-range, adatom 
diffusion) opens up. The extent to which long-range dif- 
fusion contributes to island nucleation and growth de- 
pends on the temperature. In Fig. ^, we see that at 
40 K, the DFT-kMC island density is the same as that 
for an infinitely repulsive ring (i.e., diffusing adatoms 
are unable to penetrate the ring on the time scale for 
nucleation). As the temperature increases, adatoms are 
increasingly able to penetrate the ring barrier to aggre- 
gate and add to existing islands via long-range diffusion. 
Consequently, the island density decreases with increas- 
ing temperature. It is interesting to note that for the 
conditions studied, even a relatively weak repulsive ring 
with er = 25 meV can lead to significantly higher island 
densities than those predicted by nucleation theory. 

Returning to our discussion of the experimental re- 
sults shown in Table I, we point out that the order-of- 
magnitude difference between the island densities pre- 
dicted from ab initio calculations jist and those found ex- 
perimentally for Ag on 1-ML-Ag/Pt(lll) is also seen 
in our study, comparing the island densities predicted by 
nucleation theory to those found in our DFT-kMC "com- 
puter experiments" (cf.. Fig. ||). Thus, we conclude that 
our results can explain the theoretical-experimental gap 
in the island density for Ag on 1-ML-Ag/Pt(lll). Fur- 
ther, our results indicate that for Ag on Ag(lll), the 
theoretical-experimental gap should be weaker or non- 
existent. This result is also consistent with a compari- 
son of theoretical diffusion parameters for Ag on Ag(lll) 
(vf^ = 8.2 X 10"s-\ = 82 meV) to experimen- 
tal values jTH obtained using Eq. |l| (i^o = 2 x lO^^s"^, 
= 97 meV). Finally, Bogiccvic and co-workers ||2^] 
recently found similar DFT and kMC results for both 
Al(lll) and Cu(lll) homoepitaxy. 

Thus, we conclude that long-range, electronic, 
substrate-mediated adatom interactions exist and, if 
their strength is comparable to the diffusion barrier, 
they can significantly influence surface diffusion and the 
growth morphology in thin-film epitaxy. For Ag on 
strained Ag(lll), the adatom pair interaction becomes 
repulsive past the short range and the repulsion forms a 
ring around isolated adatoms. The magnitude of the re- 
pulsion is comparable to the diffusion barrier. By inhibit- 
ing island nucleation and growth via long-range adatom 
diffusion, these interactions lead to island densities that 
are substantially larger than those predicted by nucle- 
ation theory. 
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